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In approximate Kohn-Sham density-functional theory, self-interaction manifests itself as the de- 
pendence of the energy of an orbital on its fractional occupation. This unphysical behavior translates 
into qualitative and quantitative errors that pervade many fundamental aspects of density-functional 
predictions. Here, we first examine self-interaction in terms of the discrepancy between total and 
partial electron removal energies, and then highlight the importance of imposing the generalized 
Koopmans' condition — that identifies orbital energies as opposite total electron removal energies 
— to resolve this discrepancy. In the process, we derive a correction to approximate functionals 
that, in the frozen-orbital approximation, eliminates the unphysical occupation dependence of or- 
bital energies up to the third order in the single-particle densities. This non-Koopmans correction 
brings physical meaning to single-particle energies; when applied to common local or semilocal den- 
sity functionals it provides results that are in excellent agreement with experimental data — with 
an accuracy comparable to that of GW many-body perturbation theory — while providing an ex- 
plicit total energy functional that preserves or improves on the description of established structural 
properties. 

PACS numbers: 31.15.Ew, 31.15.Ne, 31.30.-i, 71.15.-m, 72.80.Le 



I. INTRODUCTION 

Density-functional approximations) 1 " which account for 
correlated electron interactions via an explicit functional 
_E XC of the electronic density p, provide very good pre- 
dictions of total energy differences for systems with non- 
fractional occupations^ One of the notable successes of 
local and semilocal density-functional calculations is the 
accurate description of ionization processes involving the 
complete filling or entire depletion of frontier orbitals. 
In quantitative terms, the local-spin-density (LSD) ap- 
proximation and semilocal generalized-gradient approx- 
imations (GGAs) predict the energy differences of such 
reactions, namely, the electron affinity 
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(where En stands for the ground-state energy of the N- 
electron system), with a precision of a few tenths of an 
electron-volt, which compares favorably to that of more 
expensive wavefunction methods^ 

Considering the excellent performance of density- 
functional approximations in predicting total ionization 
energies, it is surprising to discover that the same theo- 
ries fail in describing partial ionization processes. As a 
matter of fact, local and semilocal functionals overesti- 
mate (by as much as 40%) the absolute energy difference 
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FIG. 1: Convexity of the LSD ground-state energy as a func- 
tion of the electron number, which results from the discrep- 
ancy between total and partial differential electron removal 
energies, i.e., X^ SD < ljf D = AjfP± < A h N S -i- 



per electron 



A 



dE 



M 



N 



dM 



(3) 



m=n+ 

of an infinitesimal electron addition and underestimate 



2 



with the same error differential energy changes 



Zjv — — 



dE 



M 



dM 



(4) 



M=N~ 



upon electron removal. 

Analytically, the discrepancy between total and partial 
differential ionization energies manifests itself into the 
convexity of E N SD as a function of N (Fig. [T]) while the 
exact ground-state energy En versus N is known to be 
described by a series of straight-line segments with pos- 
itive derivative discontinuities at integer values of N£ 
(For simplicity, we restrict the entire discussion to the 
case of the LSD functional; semilocal GGA functionals 
exhibit identical trends.) 

It is important to note that the discrepancy described 
above can be also related to the incorrect analytical be- 
havior of the LSD chemical potential /^ SD , i.e., the La- 
grange multiplier associated to particle-number conser- 
vation in the grand-canonical minimization of the total 
energy, min p {£?[/>] — p> / drp(r)}, 4-8 which can be inter- 
preted physically as the opposite electronegativity of the 
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combination of the two orthogonal ground states where 
the electron resides on either of the two ions). 

Related self-interaction artifacts explain other im- 
portant failures of local and semilocal functionals 
in predicting electron-transfer processes^ electronic 
transport^ molecular adsorption ; 14 ' 15 reaction barriers 
and energies^ - — and electrical polarization in extended 
systems^ - — Self-interaction is also connected to the 
underestimation of the differential energy gap £gap IV = 
In — An oi molecules, semiconductors, and insulators 
within LSD and GGAs^i""— 

This work is organized as follows. First, we reexamine 
the Perdew-Zunger one-electron self-interaction correc- 
tion in terms of total and differential ionization energies. 
Then, after deriving an exact measure of the unphysical 
curvature of En as a function of N based on the gen- 
eralized Koopmans' theorem, we introduce a functional 
that minimizes self-interaction errors by enforcing Koop- 
mans' condition, thereby largely eliminating the discrep- 
ancy between total and differential ionization energies 
while restoring the piecewise linearity of the ground-state 
energy En as a function of the number of electrons N. 
We conclude the study with extensive atomic and molec- 
ular calculations to demonstrate the predictive power of 
the non-Koopmans self-interaction correction. 



and determines the direction and magnitude of electron 
transfer between separated molecular fragments. 4 In fact, 
the exact dependence of the chemical potential pn when 
N varies through an integer particle number Z is known 
to be 
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where x~n denotes the Mulliken electronegativity 9 - of the 
iV-electron systemi 4 -^ 

In related physical terms, the discrepancy between ion- 
ization energies is often interpreted as arising from elec- 
tron self-interaction that causes the removal of a small 
fraction of an electron from a filled electronic state to 
be energetically less costly, in absolute value, than its 
addition to the corresponding empty state.— 

This self-interaction error is at the origin of impor- 
tant quantitative and qualitative failures that pervade 
crucial aspects of electronic-structure predictions^ Con- 
sider, for example, the dissociation of a cation dimer 
in the infinite interatomic separation limit. Here, the 
energy cost of removing a small electron fraction from 
X is lower than the energy gained in adding an electron 
fraction to X + . As a result, a portion of the electronic 
charge will transfer from X to X + , eventually leading to 
the split-charge configuration 2X2 + , with a total energy 
that is correspondingly overstabilized with respect to the 
exact solution (given, we note in passing, by any linear 



II. SELF-INTERACTION CORRECTION 

A. Perdew-Zunger one-electron correction 

Several methods have been proposed to eliminate 
self-interaction contributions and restore the internal 
consistency between total and partial electron removal 
energies. 27 ' 29 " — A widely used approach is that intro- 
duced by Perdew and Zunger^ which consists of cor- 
recting self-interaction in the one-electron approxima- 
tion by subtracting one-electron Hartree and exchange- 
correlation contributions. Explicitly, at the LSD level, 
the Perdew-Zunger (PZ) orbital-dependent functional 
and Hamiltonian are defined as 

e pz = e lsb + J2 (-E n [p i(T } - E xc [p ia ]) (7) 



= E 



LSD 



(8) 



where /9j CT (r) denotes the orbital density (a 
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represents the spin), vu(r) = SEu/8p(r) is 
the electrostatic Hartree potential, v XCy(7 (r) — 
SE xc /6p a (r) stands for the spin-dependent exchange- 
correlation potential, and the PZ one-electron 
corrective energy contributions are defined as 
Sj CT = —Enlpicr] — E xc [p ia }. 

The effect of the PZ self-interaction correction on the 
partial electron removal energies of an isolated carbon 
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FIG. 2: LSD differential ionization energies -4jv S -?i> l]f u , and 
their arithmetic mean ^(A 1 ^^ + 2^- SD ) compared with PZ 
differential ionization energies for carbon. The dotted line 
indicates the position of the experimental ionization potential. 
The behavior of a PZ functional linearly downscaled by a 
factor a pz is also shown. 



atom is depicted in Fig. [3] We observe that the LSD 
differential ionization potential 2^ SD and electron affinity 
deviate by more than 5 eV from the experimental 
total removal energy, whereas their average is in close 
agreement with experiment. The PZ correction improves 
the precision of the predicted ionization potential I^ z , 
reducing the error to less than 0.4 eV. The inaccuracy 
of the PZ electron affinity A 1 ^_ 1 remains unchanged due 
to the fact that the PZ one-electron correction vanishes 
for empty states (the slight variation of the differential 
electron affinity from 17.79 to 17.65 eV is only due to the 
self-consistent reconfiguration of the occupied manifold) . 
Consequently, the PZ functional achieves a substantial 
but only partial correction of self-interaction, reducing 
the unphysical convexity, i.e., the discrepancy An-i 
'/, -y . by a factor of \. 

Another notable feature in Fig. [5] is the overestima- 
tion of the average of the ionization potential and 
electron affinity A^L-,. Since the average ^(A^_ 1 + 
Z]y Z ) approximates J^ Z ) 10 i 40 this deviation translates into 
overestimated total ionization energies and excessively 
negative ground-state energies E^ z = —J2m=i^a?- 
(Semi-empirical relations between total energies and 
low-order ionization potentials were also evidenced by 
Pucci and March in Ref. [4l|.) In the case of iso- 
lated atoms, the energy underestimation improves (fortu- 
itously) on the description of electron correlation, bring- 
ing PZ ground-state energies in remarkable agreement 
with experiment^ However, in the more general case 
of many-electron polyatomic systems, PZ ground-state 
energies are found to be overcorrelated, yielding inac- 
curate dissociation energies and excessively short bond 
lengths 42 

Various downscaling methods have been proposed to 
correct the above trends. Nevertheless, the performance 
of such schemes is inherently limited by the fact that 
downscaling the PZ correction impairs the accuracy of 
differential ionization energies. For example, it is seen 
in Fig. [5] that a functional E aFZ , in which the self- 
interaction correction is linearly downscaled by a factor 
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5(^^-1 + 2^ Z ) without altering the preci- 
sion of 2^ z . More sophisticated downscaling methods 
exhibit identical trends i 32 i 43 ~ — 

In the next sections, we will derive a correction that 
cancels self-interaction errors for systems with fractional 
occupations without affecting the precision of total en- 
ergy differences and equilibrium structural properties for 
systems with non-fractional occupations. 



B. Measure of self-interaction 

The initial conceptual step in the construction of the 
self-interaction correction is to set forth a quantitative 
measure of self-interaction errors valid in the many- 
electron case — i.e., beyond the one-electron approx- 
imation of the Perdew-Zunger correction. With this 
self-interaction measure in hand, we will construct an 
improved self-interaction correction functional, working 
first in the simplified picture where electronic orbitals 
are kept frozen (Sec. Ill C|) . and then considering orbital 
relaxation as the next stage of refinement (Sec. Ill E[) . 

Before doing so, to derive this measure, we start from 
the physical intuition that self-interaction relates to the 
unphysical variation of the energy of an orbital ei a (/) as 
a function of its own occupation /. Hence, a necessary 
non-self-intcraction condition can be written as 
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(Here, / stands for a variable occupation while fa de- 
notes the orbital occupation that enters into the expres- 
sion of the total energy functional.) Invoking Janak's 
theorem^ this necessary condition on orbital-energy 
derivatives can be restated as a criterion on the curvature 
of the total energy: 
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where E irT (f) is the total energy E minimized under the 
constraint /j CT = / (leaving all other occupations un- 
changed). Thus, in the absence of self-interaction, the 
total energy does not display any curvature upon vary- 
ing occupations. As a corollary, any self-interaction-frec 
functional satisfies the linearity condition 



A^ ff = -f it 



dE ia (f) 



df 



0</<l, (13) 



/'=/ 



where 



AE,,, 



E ia {<d) - E la {fia) 



(14) 



4 



11,45.47 



denotes the removal energy of the orbital ipi, 

Invoking once more Janak's theorem, the above non- 
self-interaction condition yields the generalized Koop- 
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(The equivalence between Eq. (fTT]) and Eq. (fl~5j) high- 
lights the significance of Koopmans' theorem for the 
quantitative assessment of self-interaction.) It is then 
convenient to measure self-interaction in terms of the en- 
ergies II io . first introduced by Perdew and Zunger (see 
Sec. IID in Ref. [To| ) in analyzing discrepancies between 
orbital energies and vertical ionization energies^ 
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Because the self-interaction energies IIj CT quantify devia- 
tions from the Koopmans linearity [Eq. (flBf]. they will 
be termed here non-Koopmans ene rgies — the same ter- 
minology as that employed in Refs. jl0| and [39| . 
Making then use of Slater's theorem^ 
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the non-Koopmans energy can be rewritten as 
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In this form, the energy Tli a (f) is clearly seen to cor- 
respond to the integrated change of the orbital energy 
upon varying the orbital occupation — in particular, it 
is straightforward to verify that Hi a (f) = if the orbital 
energy ei a (f') does not vary with the orbital occupation 
/'• 

Using the measure defined in Eq. (fT6|) . the non-self- 
interaction criterion [Eq. (fTTj) ] can be restated exactly as 
a generalized Koopmans' condition, 
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Equation (|19[) is of central importance in this work; it 
provides a simple quantitative criterion in terms of a rig- 
orous energy-nonlinearity measure to assess and correct 
self-interaction errors. (The suggestion of using non- 
Koopmans corrections to minimize self-interaction has 
been recently introduced, in preliminary form, by Dabo, 
Cococcioni, and Marzari in Ref. pij, and, in heuristic 
form, by Lany and Zunger in Ref. 36].) The conse- 
quences of Koopmans' condition are discussed in the next 
sections. 



C. Bare non-Koopmans correction 

On the basis of the above quantitative analysis, our 
objective now is to linearize the dependence of the total 
energy as a function of orbital occupations by modifying 



the expression of the energy functional in order to cancel 
the non-Koopmans terms Hio-(f). 

To render this complex problem tractable, we first con- 
sider the restricted case where all orbitals are frozen while 
the occupation of one of them is changing in the course 
of a fictitious ionization process (the frozen-orbital ap- 
proximation). Within this paradigm, Eq. (|19[) becomes 
the restricted Koopmans ' condition, 
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where the superscript u stands for unrelaxed. Here, the 
frozen-orbital non-Koopmans energy 11^. is defined as 
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where e" CT (/) is the unrelaxed orbital energy calculated 
keeping all the orbitals frozen while setting fi a to be /. 
We underscore that in the specific case of one-electron 
systems and for / = 0, Eq. (|2"0")) yields the one-electron 
self-interaction condition of Perdew and Zunger^ 
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The restricted Koopmans' condition is thus seen to en- 
compass the Perdcw-Zungcr condition, thereby repre- 
senting a more comprehensive criterion for assessing and 
correcting self-interaction errors. 

Expectedly, satisfying the restricted Koopmans' condi- 
tion is exactly equivalent to fulfilling the restricted Koop- 
mans' theorem, 
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denotes the unrelaxed electron removal energy. 

An example of a functional that exhibits a linear 
frozen-orbital energy dependence is provided by the 
Hartree-Fock (HF) theory, generalized to fractional oc- 
cupations. Indeed, at the HF level, one can verify that 
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due to the fact that the expectation value e i 
not depend on /. 

For functionals that do not satisfy the restricted Koop- 
mans' condition [Eq. (j2"0")) ]. the unrelaxed electron re- 
moval energy can only be expressed in terms of the re- 
stricted Slater integral, 



(26) 



(note that this relation is satisfied by any functional 
whether it is subject or not to self- interaction errors). 

These considerations allow us to introduce our cor- 
rected functional aiming to satisfy Koopmans' condition; 
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at frozen orbital, it is obtained by replacing Slater terms 
- Eq. (|26p . where single-particle energies are function 
of their own occupation — with Koopmans terms - 
Eq. ([25)1 , where single-particle energies do not depend on 
occupations. In particular, we evaluate here the Koop- 
mans terms at a given orbital occupation / = / re f that 
defines the reference transition state (the determination 
of the reference occupation / ro f shall be explained on the 
basis of Slater's approximation and exchange-correlation 
hole arguments in Sec. Ill F|) . Explicitly, in the case of the 
LSD functional, the non-Koopmans (NK) self-interaction 
correction to the energy is defined as 
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where the negative sign in front of the sum follows 
from the convention that ionization energies are posi- 
tive. Rewriting Eq. (|27p in terms of the frozen-orbital 
non-Koopmans energies of Eq. (f2~Tj) , we obtain 
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where the non-Koopmans corrective terms can now be 
recast into the explicit functional form 
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where n ia {r) = |Vw| 2 (r) and p ia (r) = f ia n la (r). Here, 
the electronic density jO™ f (r) stands for the reference 
transition-state density 

pZ f (r) = f Iei n i<T (r) + £ / iff m jV (r) (30) 



P(r) + (/ref - fio)n ia (r). 



We now derive the expression of the orbital-dependent 
NK Hamiltonian. To this end, we first calculate the func- 
tional derivative of the energy contribution n"^ L D (/ re f) 
with respect to variations of the orbital density pi a . The 
expression of the functional derivatives reads 



sn. 



ti,LSD 



= (/ref - ficr)vu(r; [n i<T ]) 



Spia{r) 

+ v™»(r; [pTj]) - i£??(r; [p]) + w^%(r). (32) 
In Eq. ((32]), the potential u£ e f° CT denotes 
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where fk™^(r,r') = 8\E n + E^) / 5p a {v)5p^{v') 
is the second-order functional derivative of the LSD 
energy^ Focusing then on the cross derivatives, we ob- 
tain 
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where /^ c s ^ , cr , (r, r') is the exchange-correlation contribu- 
tion to /Hx? ffff '( r i r ')- As a nna l result, the orbital- 
dependent NK Hamiltonian can be cast into the form 



'NK 



\LSD 



-LSD 
J y:d,ia f 



(35) 



(31) where u£f ^ stands for the cross-derivative potential 
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In a nutshell, the NK Hamiltonian consists of the un- 
corrected LSD Hamiltonian calculated at the reference 
density /i LSD [p|£ f ] with the addition of two variational 
potentials. The first additional term w^fP i(T results from 
the variation of the reference density as a function of 
Pi a while the second term w]^ a springs from the cross- 
dependence of the non-Koopmans corrective terms. The 
effect of the w^^flo an d ^xdw contributions that arise 
as by-products of variationality is analyzed in the next 
section. 



D. Comparative assessment 

In this section, we assess the performance of the NK 
self-interaction correction, particularly focusing on the 
effect of variational terms on the accuracy of NK orbital 
predictions — i.e., on the cancellation of the unrelaxed 
frozen-orbital self-interaction measure n"j NK (/). 

One simple and probably the most direct way to eval- 
uate the influence of w^J^ and w]^® a is to introduce a 
non-variational orbital-energy scheme, the NKo method, 
that consists of freezing the dependence of the refer- 
ence transition-state densities and the cross-dependence 
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FIG. 3: LSD, PZ, and HF unrelaxed orbital energies and 
residual unrelaxed non-Koopmans errors of the highest oc- 
cupied state of carbon. The black arrow highlights the zero 
value for the non-Koopmans error scale. 
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puted NK and NKq orbital levels can then be compared 



for the direct assessment of w^P ia and 
plicitly, the NKq Hamiltonian can be written as 
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In the NKo optimization scheme, the Hamiltonian given 
by Eq. (f38|) is employed to propagate orbital degrees of 
freedom at fixed p™ . Reference transition-state densities 
are then updated according to Eq. (|30l) . The procedure 
is iterated until self-consistency. 

Due to the loss of variationality, the obvious practi- 
cal limitation of the non-variational NKo orbital-energy 
method is that it cannot provide total energies and inter- 
atomic forces. However, NKo is of great utility in evalu- 
ating the intrinsic performance of the NK correction. In 
itself, the NKo formulation is also useful in determining 
orbital energy properties that are particularly affected by 
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Focusing now on computational predictions, the occu- 
pation dependencies of the LSD, HF, PZ, and NK unre- 
laxed orbital energies 
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FIG. 4: Non-Koopmans unrelaxed orbital energies and resid- 
ual unrelaxed non-Koopmans errors of the highest occupied 
state of carbon for select arbitrary values of the reference oc- 
cupation (/ rc f = i, |, and |) using (a) the NK total-energy 
method and (b) the NKo orbital-energy method. The black 
arrow highlights the zero value for the non-Koopmans error 
scale. 



of corrective energy terms, thereby eliminating i(T 
and ?ct contributions to the effective potential. Com- 



of the highest atomic orbital of carbon are depicted in 
Figs. [3] and |3U a). The salient feature of the LSD graph is 
the large variation of the orbital energy from —19.40 to 
—6.15 eV, reflecting the strong nonlinearity of the corre- 
sponding unrelaxed ionization curve. The PZ variation 
is found to be twice lower than for LSD, confirming the 
trends observed in Sec. Ill Al In contrast, the HF unre- 
laxed ionization curve exhibits a perfectly linear behav- 
ior (i.e., the unrelaxed orbital energy remains constant). 
This trend is closely reproduced by the NK functional 
[Fig. Ufa)]; on the scale of LSD residual non-Koopmans 
errors, the eye can barely distinguish any deviation of 
the NK unrelaxed orbital energies as a function of fi a 
regardless of the value of the reference occupation. 

The above observation is due to the fact that the vari- 
ational contribution w^f^ affects the orbital energy efj 1 
indirectly, i.e., only through the self-consistent response 
of the orbital densities since {^i^w^^^ia) = at self- 
consistency. Furthermore, a Taylor series expansion of 
*xd,?o- rev eals that (^l^xf ^Jfe) does not cause no- 
table departure from the linear Koopmans behavior. In 
quantitative terms, the dominant term in the expansion 
of the residual NK self-interaction error is of the fourth 
order in orbital densities: 
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(where /xc,cti 2 ...„ ( r i2 - n) denotes the nth order func- 
tional derivative of the LSD exchange-correlation en- 
ergy), whereas the PZ correction is found to be less ac- 



Despite the very good accuracy of the NK correction, 
direct confrontation with NKo results — for which the 
non-Koopmans measure 

Ka (/) [Eq- (EH)] is obviously 
canceled for any value of / — reveals that NK tends 
to underestimate orbital energies with deviations of 0.5 
to 1.5 eV that gradually increase with / re f [Fig. HJa,b)]. 
The practical consequences of this observation will be 
discussed in Sec. IIII1 

As a conclusion of this preliminary performance eval- 
uation, the NK frozen-orbital correction results in a con- 
siderable reduction of residual errors IL^. (/), bringing 
density-functional approximations in nearly exact agree- 
ment with the frozen-orbital linear trend while exhibit- 
ing a slight tendency to underestimate orbital energies. 
In the next sections, we present an extension of the NK 
correction beyond the frozen-orbital paradigm. 

E. Screened non-Koopmans correction 

Having derived and examined the bare non-Koopmans 
frozen-orbital correction, we now turn to the analysis of 
self-interaction for relaxed orbitals. The effect of the NK 
correction on relaxed partial ionization energies for car- 
bon is illustrated in Fig. [SJ The first important obser- 
vation is that the NK correction decreases the unphysi- 
cal convexity of En, reducing the unphysical discrepancy 
(An-i — %n) by a factor of 6 regardless of the reference 
occupation / rc f. 

The second notable feature is the fact that the NK cor- 
rection reverts the convexity trend (i.e., 2 N K > A^fz-y), 
transforming the convex dependence of En into a piece- 
wise concave curve. In fact, bearing in mind that relax- 
ation contributions to En are always negative, it is rel- 
atively straightforward to show that any functional that 
satisfies the restricted Koopmans' condition [Eq. (|20p] is 
piecewise concave. Consequently, in contrast to the aPZ 



(40) 

I 

curate in minimizing the self-interaction measure by one 
order of precision: 



(41) 

I 




eV 



eV 



FIG. 5: LSD differential ionization energies .4jv_i, X N , and 
their average ^(A]f^i +Ijv SD ) compared with NK differential 
ionization energies for (a) / re f = | and (b) / re f = | for carbon. 



functional introduced in Sec. Ill A| there always exists a 
value of the coefficient a NK (0 < a NK < 1) for which the 
oNK functional 

E aNK = E LSU + a NK (£ NK - £ LSD ) (42) 

restores the agreement between A N K _ 1 and T N K - At this 
crossing point, the aNK relaxed ionization curve is in 
close agreement with the exact linear trend described by 
the generalized Koopmans' condition (see Sec. IIII A[) . 
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Making the approximation that ionization energies 
vary linearly with a NK , the value of the coefficient for 



which the crossing A% N J^ 
mated as 



1 N NK occurs can be esti- 



,NK 



4 LSD _ tLSD 



(AhSD _ TLSD>> _ (A 
V^N—l N I 



NK 



tNK 

N 



(43) 



This initial estimate can then be refined using the secant- 
method recursion 



(1 



Otn+l 



\i io„NK 



■ 1 



■a„NK 
'jV 



( <o„NK 



Tti„NK 



) 



- (A NK 



T NK\ 



(44) 



In practice, it is observed that only one or two iterations 
are sufficient to determine the value of the coefficient 
q nk _ ii rtln _ >oo Q, n anc [ bring the difference A^J^ —X N NK 
below 0.2 eV. 

Physically, the coefficient a NK is directly related to the 
magnitude of orbital relaxation upon electron removal. 
Indeed, a NK can be viewed as a screening coefficient 
whose value is close to I for weakly relaxing ionized sys- 
tems (since E^ K is already piecewise linear in the absence 
of orbital relaxation) and is found to be small when re- 
laxation is strong (see Sec. IIII A[) . The corrective factor 
a NK introduced here is thus endowed with a clear phys- 
ical interpretation. 

As a final note, it should be emphasized that we 
have adopted here a simple picture of orbital relaxation 
through the coefficient a NK , which can be viewed as 
a uniform and isotropic screening factor. More elabo- 
rate screening functions could be employed (at the price 
of computational complexity). Such accurate screen- 
ing approaches provide promising extensions of the aNK 
method and represent an interesting subject for future 
studies. 



into the definition of the non-Koopmans energy contri- 
butions [Eq. (|2ip], one can demonstrate that 



nU,LSD/ r 
N Uref 



to arrive at the relation 



t-qNK 
1 N 



(/re 



d 



j)«0 



rLSD 
N ■ 



(48) 



(49) 



This result explains the accuracy of the aNK correction 
in predicting total electron removal energies with the ref- 
erence occupation / re f = \. 

It should be noted that in the particular case of one- 
electron systems, / re f = | is not the only possible choice 
for the reference occupation. Indeed, due to the fact the 
Hartree and exchange-correlation potentials vanish for 
empty solitary orbitals, the approximation 



AE 



LSD 



m,LSD 



(/ = o) 



(50) 



also holds; hence, the solution / re f = is equally valid. In 
fact, / re f = leads to the exact one-electron Hamiltonian 
and the exact solution to the one-electron Schrodinger 
problem. 

As an alternative to transition-state arguments, the 
value of / rc f can be justified by inspecting the sum 
rule satisfied by the exchange-correlation hole (xc- 
hole) j 10 ' 50 ' 51 The xc-hole h xc is defined by the relation 



1 f drdr> P(r)h xc (rr>) 



(51) 



and can be explicitly written through the adiabatic con- 
nection formalism as 

p(r)/ lxc (r,r')-(/5 2 (r,r')) A -p(r)p(r') (52) 



F. Reference transition-state occupation 



with 



We now proceed to examine the influence of the ref- 
erence occupation on the accuracy of calculated total 
electron removal energies. To this end, we compare the 
average of the aNK differential electron removal ener- 
gies (that closely approximates the total ionization po- 
tential) of carbon for / re f different from and equal to | 
in Figs. [UJa) and[3]Jb), respectively. In the former case, 
the diagram indicates that the average 



1 



/ aoNK , 7", 



N 



mNK 
2 N 



(45) 



deviates significantly from its LSD counterpart, whereas 

l 

2 ' 



such deviations do not occur in the case / rc f = -k- In fact, 



from Eqs. ©, ([25]). and ([32"f. it can be shown that 

raNKfr \ ^ tLSD NKttM.LSD/,. n , 

(neglecting orbital reconfiguration). Then, substituting 
the (restricted) Slater's approximation^ 

LSD u,LSD 



AE, 



{ficr 



(47) 



(p 2 (r,v')) x = / dA(* A |^ 2 (r,r')|*A>. (53) 



In Eq. ([521), = ^t( r )^t( r ')^,( r ')^( r ) denotes the 

pair density operator and ^ \ stands for the ground state 
of a fictitious system where the Coulomb interaction is 
scaled down by A and the effective local single-electron 
potential v\ is added to the Hamiltonian to keep the 
ground-state density constant with respect to A. 

By definition , 10 ! 52 the exact xc-hole corresponding to 
a system with a fractional number of electrons N + uj, 
satisfies the sum rule 



f dr'h xc (r,r') = -1 



+ 0,(1-0,) / V^+iW-Pa.atW (54) 
Jo P(r) 

where the ground-state density p\ t m of the system with 
M = N and M = N + 1 electrons depends on A due 
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to the fact that M differs from N + oj. It is important 
to note that the exact xc-hole h xc integrates to —1 only 
for integer electron numbers at variance with the LSD 
xc-hole /ixf D , which integrates to —1 irrespective of the 
number of electrons. 

Turning now to the aNK xc-hole h^ K , it is possible 
to derive a relation similar to Eq. (|54|) : 

J dr'^ c NK (r,r') = -1 

+ a NK £iW( 2 ^-/-) ! ¥r (55) 

^ P r 



for N > 1 (the detailed derivation is presented in Ap- 
pendix EJ. It is then clear that the occupation / rc f = i 
allows to satisfy the xc-hole sum rule exactly for integer 
number of electrons, and at least approximately for frac- 
tional numbers. (Note that in the case N < 1, the value 
/rcf = corresponds also to the exact sum-rule.) For 
PZ, it has been argued that the enforcement of a similar 
sum rule is critical to the quality of the self-interaction 
correction. 4s 



As a result, the explicit expression of the screened aNK 
functional with / re f = ^ reads 



E 



aNK 



E 



LSD 



+ « NK E 



fi,(l-fia)Ea[r. 



drv^V I r: 



P+ _ fi<r)lH<J 



p ia {v)-E^[p] 



(56) 



In the next section, we assess the predictive accuracy of 
the aNK functional given by Eq. (|55|) . first focusing on 
atoms then extending applications to molecular systems. 



III. APPLICATIONS 
A. Atomic ionization 

In order to probe the performance of the aNK func- 
tional, we calculate the electron removal energies of a 
complete range of atomic elements, from hydrogen to 
xenon, using the all-electron LD1 code of the Quantum- 
Espresso distribution^ (It should be noted that the re- 
ported atomic calculations do not account for relativistic 
effects; the inclusion of scalar-relativistic contributions 
results in marginal deviations of <80 meV in the abso- 
lute precision of aNK differential electron removal en- 
ergies.) The LD1 code proceeds by iterative integration 
of the spherically symmetric electronic-structure prob- 
lem on logarithmic grids. The validity of this integra- 
tion procedure, also employed by Perdew and Zunger in 
their original work,— has been extensively discussed and 
carefully verified by Goedecker and Umrigar. 42 More re- 
cently, Stengel and Spaldin have shown that the spher- 
ical approximation is actually required to prevent the 
appearance of unphysical hybrid states in the electronic 
structure of isolated atoms and predict physical atomic 
electron removal energies within orbital-dependent self- 
interaction corrections^ 

We first compute the screening coefficient a NK 
[Eqs. (H3J) and (05)] for the elements of the five first pe- 
riodic rows based on the electronic configurations tabu- 
lated in Ref. [54j (with the exception of Ni for which we 



used the slightly unfavored 4s 1 3d 9 configuration instead 
of 4s 2 3d 8 that overestimates the experimental ionization 
potential by more than 1 eV) . The calculated a NK are re- 
ported in Fig. [6] The graph confirms that the screening 
coefficient varies in a narrow range. Not unexpectedly, 
the value of a NK is found to be maximal for the elements 
whose outermost electronic shell contains a single elec- 
tron, such as hydrogen for which a NK = 1 and alkali 
metals, e.g., a NK (Na) = 0.99. In contrast, the coeffi- 
cient a NK saturates to low values for filled shell elements, 
namely, noble gases, e.g., a NK (Hc) = 0.66, alkaline earth 
metals, e.g., a NK (Be) — 0.72, and filled shell transition 
metals, e.g., a NK (Pd) = 0.69. 

After calculating atomic screening coefficients, we com- 
pare aNK differential electron affinity predictions with 
LSD and experiment in Fig.|S^ The comparison demon- 
strates the predictive ability of the aNK method, which 
brings partial electron removal energies An-i in very 
close agreement with experimental total electron removal 
energies A c ^ 1 , whereas LSD is found to considerably 
overestimate An~i- In quantitative terms, the differen- 
tial LSD energy A]^^ is overestimated by more than 
4 eV with a standard deviation of 1.85 eV [Fig. EJa)]. 
Comparable deviations are obtained with the PZ self- 
interaction correction. The HF energy A^f_ 1 are instead 
underestimated by a smaller margin of 1.48 eV. The aNK 
correction results in substantial improvement in the cal- 
culation of partial electron removal energies, reducing the 
error to 0.31 eV. Here, it is quite interesting to note 
that the aNK variational contributions counterbalance 
the slight tendency of the aNKo correction to underesti- 
mate electron removal energies within LSD. 

We now examine partial ionization potential predic- 
tions (Fig. O. A marked difference with the above elec- 
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FIG. 6: Non-Koopmans screening coefficient a NK for the eiements of the five first periods 
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FIG. 7: Atoms: mean deviation (MD), mean absolute devia- 
tion (MAD), and root mean squared deviation (RMS) of LSD, 
HF, PZ, aNK, and aNKo atomic differential electron removal 
energies (a) An-i and (b) In relative to experiment for the 
elements of the five first periods. Energy errors are in eV. 



tron affinity results is the enhanced accuracy of the HF 
and PZ theories. The improved performance in predict- 
ing atomic ionization potentials results from the fact that 
orbital relaxation compensates the absence of correlation 



contributions in HF and cancels residual non-Koopmans 
errors in PZ^ Nevertheless, even with beneficial error 
cancellation in favor of HF and PZ, the aNK deviation 
is still the lowest, approximately equal to the PZ mean 
absolute error of 0.34 eV [Fig. [7(b)]. 

The precision of aNK differential ionization energies 
reflects the intrinsic accuracy of the underlying LSD to- 
tal energy functional in reproducing subtle atomic ion- 
ization trends that are difficult to describe with, e.g., 
semi-empirical approximations^ 

To complement these observations, the full dependence 
of the LSD, PZ, and aNK relaxed energies £2 P of the 
highest occupied atomic orbital of carbon as a function 
of its occupation is depicted in Fig. [TUl Confronting the 
LSD and PZ graphs with those presented in Fig. [3[ it is 
seen that orbital relaxation causes a non-negligible de- 
crease of the unphysical shift £2 P (1) — £2 P (0) of 1.5 eV 
for both functionals. Additionally, the PZ orbital energy 
e 2p becomes less curved at higher occupations, confirm- 
ing that orbital relaxation enhances the performance of 
the PZ correction^ Nevertheless, the inflexion of the 
curve remains important in the vicinity of /2 P = due to 
the fact that the PZ correction leaves the energy of the 
empty state unchanged. We observe that this unphysical 
trend is almost completely removed by the aNK correc- 
tion, clearly showing that the screened aNK method is 
apt at imposing the generalized Koopmans' condition for 
any fractional value of /2 P . 

The above comparisons demonstrate the predictive 
performance of the non-Koopmans method in correcting 
atomic differential electron affinities and first ionization 
potentials, placing An-i and In predictions on the same 
level of accuracy with respect to experiment. The fact 
that aNK improves An-i and Zjy with the same preci- 
sion ensures the accuracy of non-Koopmans total energy 
differences and related equilibrium properties. The re- 
sults presented in the next sections provide further sup- 
port to this conclusion. 
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FIG. 8: LSD and «NK differential electron affinities An-i (i.e., opposite energy of the lowest unoccupied orbital of the ionized 
atom X + ) compared with experiment for the elements of the five first periods. 
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FIG. 9: LSD and aNK differential ionization potentials Xn (i.e., opposite energy of the highest occupied orbital of the neutral 
atom X) compared with experiment for the elements of the five first periods. 
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FIG. 10: LSD, PZ, and aNK relaxed orbital energies of the 
highest occupied state of carbon. 
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FIG. 11: Molecular non-Koopmans screening coefficient a 



B. Molecular ionization 

In this section, we focus on the study of molecular sys- 
tems. For this purpose, we have implemented the HF, 
PZ, and aNK methods in the plane-wave pseudopoten- 
tial CP (Car-Parrinello) code of the Quantum-Espresso 



distribution. 53 In this code, orbital optimization proceeds 
via fictitious Newtonian damped electronic dynamics. 

The main difficulty in the CP implementation of the 
HF, PZ, and aNK functionals is the correction of 
periodic-image errors that arise from the use of the su- 
percell approximation^ Such numerical errors preclude 
the accurate evaluation of exchange terms and orbital 
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TABLE I: LSD, HF, PZ, aNK, and qNKo differential molecular electron removal energies compared with experimental vertical 
electron removal energies. Mean deviations (MD), mean absolute deviations (MAD), and root mean squared deviations of 
the error (RMS) in absolute and relative terms are also reported. The adiabatic ionization lower bound is given when the 
experimental vertical ionization energy is not available. Energies are in eV. 



LSD 


HF 


PZ 


aNK 


qNK 


Exp. a 


An-i In 


An-i In 


An-i In 


An-i In 


An-i In 





H 2 


18.84 


10.16 


N 2 


20.85 


10.37 


2 


18.94 


7.20 


P2 


14.28 


7.26 


S2 


13.28 


5.81 


PH 


13.92 


5.81 


HC1 


17.81 


8.11 


CO 


18.70 


9.14 


cs 


15.20 


7.42 


H 2 


18.97 


7.33 


H 2 S 


14.87 


6.39 


NH 3 


16.07 


6.23 


PH 3 


14.50 


6.84 


CH4 


18.69 


9.46 


SiH 4 


15.86 


8.50 


C2H2 


16.19 


7.40 


C2H4 


15.12 


7.01 


MD 


4.57 


-4.35 




38.5% 


-36.4% 


MAD 


4.57 


4.35 




38.5% 


36.4% 


RMS 


0.95 


0.63 




7.7% 


4.0% 



14.78 


16.24 


19.01 


12.77 


17.15 


21.55 


9.27 


14.71 


18.11 


8.56 


10.84 


14.40 


7.49 


10.49 


13.05 


8.39 


10.25 


13.81 


10.33 


13.05 


17.58 


11.03 


15.06 


18.45 


7.30 


12.83 


14.53 


8.97 


13.81 


18.50 


8.17 


10.55 


14.70 


7.61 


11.55 


15.75 


8.29 


10.46 


14.15 


11.95 


14.92 


18.98 


10.84 


13.26 


16.46 


8.79 


11.40 


16.19 


7.91 


10.43 


15.07 


-2.46 


0.75 


4.47 


-21.0% 


5.9% 


37.4% 


2.46 


0.80 


4.47 


21.0% 


6.4% 


37.4% 


0.83 


0.73 


0.89 


7.4% 


5.9% 


6.3% 
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10.17 


10.19 


10.84 
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10.02 


13.82 


13.09 


13.17 


15.40 


14.15 


14.24 


13.29 


11.63 


11.76 


14.77 


13.17 


13.45 


11.52 


10.77 


10.86 


12.50 


11.15 


11.40 


11.45 


10.77 


10.83 


16.20 


14.46 


14.51 


14.35 


12.67 


12.68 


12.97 


12.05 


12.14 
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11.43 
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0.40 


0.49 


13.7% 


3.4% 


4.2% 


1.66 


0.50 


0.58 


13.7% 


4.1% 


4.8% 


0.66 


0.46 


0.48 


5.0% 


3.7% 


3.9% 



15.74 


15.69 


>15.43 


15.38 


15.52 
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12.30 


10.35 


10.35 


10.62 


9.59 


9.54 


9.55 


9.68 


9.68 
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12.30 


12.32 
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13.68 


13.74 


14.01 


10.83 


10.89 


>11.33 


11.75 


11.94 
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10.15 


10.17 
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10.18 


10.31 
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10.41 


10.41 


10.59 


13.91 
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11.49 
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0.40 
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electrostatic potentials. To eliminate periodic-image er- 
rors in the plane- wave evaluation of exchange and elec- 
trostatic two-electron integrals, we employ countercharge 
correction techniques^ In addition to this difficulty, ex- 
plicit orthogonality constraints must be considered for 
the accurate calculation of the gradient of the orbital- 
dependent PZ and aNK functionals42 To incorporate 
these additional constraints, we use the efficient itera- 
tive orthogonalization cycle implemented in the original 
CP code. 57 In terms of computational performance, the 
cost of aNK calculations is here only 40% higher than 
that of PZ and lower than that of HF. 

In Table [fl we compare LSD, HF, PZ, and NK par- 
tial electron removal energy predictions for a representa- 
tive set of molecules. In each case, molecular geometries 
are fully relaxed (the accuracy of equilibrium geometry 
predictions will be examined in Sec. 1111 CI) . To per- 
form our calculations, we employ LSD norm-conserving 
pseudopotentials 5 -^ with an energy cutoff of 60 Ry for 
the plane- wave expansion of the electronic wavefunctions. 
With this calculation parameter, we verify that *4/v-i 
and In are converged to within less than 50 meV. 



It is frequently argued that substituting LSD pseu- 
dopotentials for their HF, PZ, and NK counterparts has 
minor effect on the predicted energy differences. 42 Com- 
paring our pseudopotential calculations with all-electron 
atomic results (see Sec. IIII Al) . we actually found that 
the use of LSD pseudopotentials yields HF, PZ, and NK 
electron removal energies with a typical error of 0.1 to 0.2 
eV. However, since these moderate deviations affect HF, 
PZ, and NK predictions in identical manner, the pseu- 
dopotential substitution does not alter the validity of the 
present comparative analysis. 

As expected, one conspicuous feature in Table U is the 
poor performance of LSD that predicts molecular partial 
electron removal energies with an average error of ±40%. 
As was the case for atoms, the PZ self-interaction correc- 
tion reduces the error in predicting Ijv to less than 14%, 
which corresponds to an average deviation of 1.68 eV, 
whereas An-i predictions are not improved. In com- 
parison, aNK partial ionization energies are predicted 
with a remarkable precision of 0.50 eV (4.1%) and 0.58 
eV (4.8%) for A N -i and T N , respectively. The aNK 
accuracy in predicting molecular vertical ionization ener- 
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gies compares favorably (arguably, even more accurately) 
with that of recently published fully self-consistent GW 
many-body perturbation theory calculations. 59 

It should also be noted that we perform here only one 
iteration [Eq. (|43l) ] to determine the screening coefficient 
and that the calculated a NK vary in a very limited range 
of values, even narrower than that found in the case of 
atoms (Fig. [TTj) . 

To conclude the analysis of An-i and In predictions, 
we focus on the influence of the variational contributions 
and Wxd?cr on * ne accuracy of the aNK total en- 



.-.LSD 

U Tcf,i<7 



ergy method. Similarly to the comparative analysis pre- 
sented in Sec. Ill Dt we use the aNKo orbital-energy for- 
mulation to evaluate the magnitude of w^J 1 ^ and 
errors. In agreement with the trend already observed, 
aNKo predictions for An-i and In are lower than aNK 
electron removal energies with negative shifts of 0.59 eV 
and 0.68 eV, respectively. (Thus, aNKo results are found 
to be even closer to experiment than their aNK coun- 
terparts with mean absolute error margins of only 3.2% 
for An-i and 2.5% for In-) This direct comparison in- 
dicates that u>ref?CT an d ^xd?CT introduce non-negligible 
energy shifts in the calculation of frontier orbital lev- 
els. However, these errors are much smaller than typical 
self-interaction deviations of 4 to 5 eV, providing quanti- 
tative justifications of the excellent precision of the aNK 
method. 



C. Equilibrium structural properties 
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FIG. 12: LSD, PZ, and aNK molecular bond lengths com- 
pared with experiment (Ref. [6(|). 

After analyzing partial electron removal energies, we 
now report on the accuracy of aNK equilibrium geometry 
calculations. We compare LSD, PZ, and aNK structural 
predictions to experimental bond lengths in Fig. [T2] and 
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FIG. 13: Molecular geometries: LSD, HF, PZ, and aNK mean 
deviations (MD), mean absolute deviations (MAD), and root 
mean squared deviations of the error (RMS) with respect to 
experiment of the predicted bond lengths for the 17 molecules 
listed in Tabled] Error bars are in A. Relative errors are also 
reported. 



we present LSD, HF, PZ, and aNK error bars in Fig. IT3"1 
The first important observation is the very good accuracy 
of LSD predictions with a mean absolute relative error 
of 1.1% for the seventeen molecules listed in Table HI PZ 
bond lengths are instead sensibly underestimated with a 
mean uncertainty of 2.8%. In contrast with PZ calcula- 
tions, aNK results deviate from experiment by a relative 
error margin of 0.8%, which is lower than that of LSD, 
demonstrating that the aNK self-interaction correction 
does not deteriorate and even improves LSD structural 
predictions, at variance with the conventional PZ self- 
interaction correction. 

These results illustrate the tendency of PZ to overbind 
molecular structures, and confirm the systematic im- 
provement brought about by the aNK correction. The 
promising potential of the aNK correction in predicting 
other thermodynamical properties (e.g., dissociation en- 
ergies and vibrational frequencies) will be critically ex- 
plored in a separate study. 



D. Photoemission energies 

Having validated the non-Koopmans self-interaction 
correction for the calculation of electron removal ener- 
gies and equilibrium structures, we now evaluate the per- 
formance of the aNK and aNKo methods in predicting 
photoemission energies, for which LSD and GGAs exhibit 
notable failures. 

From the theoretical point of view, the very poor per- 
formance of LSD and GGA is expected; Kohn-Sham 
density-functional theory eigenvalues are not meant to 
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TABLE II: LSD, HF, PZ, and aNK orbital energies of neon, 
argon, and krypton compared with experimental photoemis- 
sion energies. Relative mean absolute deviations (MAD) with 
respect to experimental photoionization results are also re- 
ported. The experimental photoemission energies of the spin- 
orbit doublets of p and d orbitals are indicated. Computa- 
tional photoionization predictions do not include spin-orbit 
coupling. Energies are in eV. 





LSD 


HF 


PZ 


aNK 


Exp. a 


Ne 2p 


13.54 


23.11 


22.91 


22.52 


21.6-21.7 


2s 


35.99 


52.49 


45.13 


45.11 


48.5 


Is 


824.68 


891.75 


889.41 


872.14 


870.2 


I\T op 


i n /in 


10. uo 


10. i D 


1D.U4 


1 C 7 ICQ 

io. / — io. y 


3s 


24.03 


34.74 


30.22 


30.54 


29.3 


2p 


229.77 


260.45 


256.12 


254.65 


248.4-250.6 


2s 


293.73 


335.30 


315.49 


315.40 


326.3 


Is 


3096.69 


3227.47 


3218.88 


3193.55 


3205.9 


Kr 4p 


9.43 


14.25 


13.97 


14.35 


14.1-14.2 


4s 


22.33 


31.34 


27.78 


28.27 


27.5 


3d 


83.65 


104.06 


101.29 


101.67 


93.8-95.0 


3p 


192.84 


226.70 


209.04 


210.71 


214.4-222.2 


3s 


253.48 


295.21 


269.48 


271.24 


292.8 


2p 


1633.17 


1714.53 


1695.09 


1692.46 


1678.4-1730.9 


2s 


1803.75 


1902.11 


1852.00 


1853.32 


1921.0 


Is 


13877.37 


14154.31 


14128.17 


14080.07 


14326.0 


MAD 


19.2% 


4.5% 


3.3% 


3.2% 
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TABLE III: LSD, HF, PZ, and aNK orbital energies of 
benzene compared with experimental photoemission energies. 
Relative mean absolute deviations (MAD) with respect to ex- 
perimental photoionization results are also reported. Energies 
are in eV. 



LSD HF PZ aNK aNK Exp. 



ei g 


6.59 


9.18 


9.43 


10.39 


9.39 


9.3 


e 2g 


8.28 


13.54 


15.46 


12.66 


12.48 


11.8 


a2u 


9.43 


13.64 


12.99 


13.25 


12.60 


12.5 


eiu 


10.33 


16.02 


17.67 


14.75 


14.56 


14.0 




11.02 


16.95 


18.40 


15.46 


15.15 


14.9 


bi u 


11.26 


17.51 


18.82 


15.65 


15.69 


15.5 


ai 9 


13.10 


19.26 


20.60 


17.58 


17.30 


17.0 


e 2g 


14.85 


22.39 


22.72 


19.27 


19.34 


19.2 


MAD 


26.1% 


12.0% 


18.4% 


4.9% 


2.1% 
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predict excited-state properties j ' 65 In practice, total 
electron removal energies computed from constrained 
density-functional calculations (A C SCF) are typically 
found to be in good agreement with experiment j 62 i 66 
This level of accuracy suggests that orbital-dependent 
self-interaction corrected functionals can provide orbital 
energies in accordance with spectroscopic results^ This 
expectation is confirmed by PZ photoemission predic- 
tions for neon, argon, and krypton that we reproduce 
here using the LD1 code (Table HT1). 

Nevertheless, similarly to the trend observed in Sec. 
IIIIB1 the predictive ability of PZ deteriorates in the case 
of molecular photoionization; this is at variance with the 
aNK method. To illustrate this fact, we compare LSD, 
HF, PZ, and non-Koopmans predictions for the photoe- 
mission spectrum (PES) of benzene in Table IIIII and 
fullerene in Table IIV1 In the molecular photoemission 
calculations, we use the CP code with the computational 
procedure described in Sec. IIIIBI We employ fully re- 
laxed geometries for benzene and the LSD atomic struc- 
ture of Ceo, which is found to be in excellent agreement 
with the NMR experimental geometry^ 

Focusing first on benzene, we observe that LSD under- 
estimates electron binding energies with errors as large 
as 4.35 eV for low-lying states. In contrast to LSD, the 
HF theory provides overestimated photoemission ener- 
gies with absolute deviations that increase gradually from 
0.12 to 3.19 eV when approaching the bottom of the PES. 
Similar trends are observed for PZ with the difference 
that the errors do not systematically increase with in- 
creasing photoemission energies, leading in particular to 
the incorrect ordering of the ei g and a,i u levels. In con- 
trast, aNK restores the correct relative peak positions 
and yields slightly overestimated electron binding ener- 
gies with an absolute precision of 4.9%. The slight ten- 
dency of aNK to overestimate electron binding energies 
is here again due to the influence of variational contri- 
butions, as directly confirmed by the performance of the 
aNKo orbital-energy method, which predicts photoemis- 
sion energies in remarkable agreement with experiment. 

Similarly to benzene, LSD energy predictions for 
fullerene are significantly underestimated. However, 
since the dispersion of the errors is much narrower than 
in the case of benzene, a simple shift of LSD photoe- 
mission bands, equal to the difference between the the- 
oretical and experimental HOMO levels, can bring the 
predicted PES in close agreement with experiment^ De- 
spite the excellent precision of HF in the top region of the 
spectrum, HF photoemission energies are largely overes- 
timated for low-lying states. In addition, HF inverts the 
h g and g g states in the second photoemission band al- 
though it predicts the correct peak ordering in the third 
and fourth bands f& The performance of PZ is found to 
be slightly worse than that of HF with significant qualita- 
tive errors in the grouping and ordering of the states. In 
contrast, aNK correctly shifts the spectrum and brings 
photoemission energies in very good agreement with ex- 
periment. Predicted aNK binding energies are also in 
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TABLE IV: LSD, HF, PZ, and aNK orbital energies of fullerene Ceo compared with constrained LSD total energy differences 
(A C LSD) and experimental photoemission energy bands. Energies are in eV. 



band 




LSD 




HF 




PZ 




aNK 




A c LSD a 


Exp. 6 


I 


llu 


5.84 




7.49 


llu 


8.77 


llu 


7.45 


llu 


7.61 


7.60 


II 


gg 


7.03 




9.42 


gg 


9.80 


g 9 


8.64 


g 9 


8.78 


8.95 






7.15 


gg 


9.64 


li S 


10.48 


h 9 


8.75 


h 9 


8.90 




III 


K 


8.72 


gu 


12.42 


gu 


12.21 


llu 


10.31 


llu 


10.47 


10.82-11.59 




gu 


8.74 


tj U 


12.99 


t U 


12.69 


gu 


10.35 


gu 


10.50 








9.03 


llu 


13.08 






h 9 


10.64 


h 9 


10.79 








9.28 










t u 


10.91 




11.03 




IV 


gu 


10.05 


K 


13.46 


llu 


14.13 


g« 


11.66 


gu 


11.79 


12.43-13.82 




t 9 


10.52 


gu 


15.06 


t 9 


15.41 


t 9 


12.12 


t 9 


12.28 






h 9 


10.59 


t 9 


15.20 


h 9 


15.81 


h 9 


12.20 


h 9 


12.33 










h 9 


15.66 
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been demonstrated for a range of atomic and molecular 
systems. 

The theory developed here represents a significant step 
in the correction of electron self-interaction in electronic- 
structure theories. Nevertheless, interesting problems are 
left open. One central question is that similarly to the PZ 
approach, the aNK method leads to an orbital-dependent 
Hamiltonian, although it is always possible in principle to 
derive a consistent density-dependent formulation using, 
e.g., optimized effective potential mappings<21 It is a long 
held tenet that the orbital dependence of self-interaction 
functionals and the subsequent loss of invariance with re- 
spect to unitary transformation of the one-body density 
matrix precludes applications to periodic systems (e.g., 
conjugate polymers and crystalline materials). In fu- 
ture studies, we will explore solutions to this central con- 
ceptual difficulty without resorting to density-dependent 
unitary invariant mappings. 
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excellent agreement with constrained LSD total energy 
differences f2S providing a final validation of the perfor- 
mance of the aNK self-interaction correction in bringing 
physical meaning to orbital energies — i.e., in identifying 
orbital energies as opposite total electron removal ener- 
gies. 



IV. SUMMARY AND OUTLOOK 

In summary, we have demonstrated that the correc- 
tion of the nonlinearity of the ground-state energy En 
as a function of the number of electrons N, at the origin 
of important discrepancies between total and differential 
electron removal energies, and related fundamental qual- 
itative and quantitative self-interaction errors, can be 
achieved without altering the otherwise excellent perfor- 
mance of density-functional approximations in describ- 
ing systems with non-fractional occupations. To con- 
struct the non-Koopmans self-interaction correction, we 
have first defined an exact non-Koopmans measure of 
self-interaction and adopted the frozen-orbital approxi- 
mation (i.e., the framework of the restricted Koopmans' 
theorem) as a working alternative to the conventional 
one-electron paradigm. We have then accounted for or- 
bital relaxation by introducing the screening coefficient 
a NK , which bears the physical significance of a uniform 
and isotropic screening factor that can be determined 
itcratively — thereby closely satisfying the generalized 
Koopmans' condition. This self-interaction correction 
scheme can be applied to any local, semilocal or hy- 
brid density-functional approximation. The remarkable 
predictive performance of the non-Koopmans theory has 
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Appendix A: Screened non-Koopmans exchange-correlation hole sum rule 

In this appendix, we derive the explicit expression of the aNK xc-hole. Starting from the relation 

E aNK = £ LSD + « NK £ (i» - pu,] - I»] + f drp ia (r)v^ <<T (r; \f£\)) (Al) 



and from the definition of the xc-hole [Eq. (|5ip], the contributions to the total aNK xc-hole arising from the three 

xc.icr' xc.zcr' 



summation terms in Eq. (|A1|) can be worked out. Those terms will be labeled h^] ia , h^ ia , and , respectively. 



Focusing on the first term, it is straightforward to obtain 



C( r > r ) 



p(r) ' " p(r) 



-(p(r')-p ia (r')) + 



h^(v,v';[p-p ia ]). 



(A2) 



Turning to the second term and including the appropriate sign, one obtains 



4 b c ^(r, r') = -p(r') - h^(r, r'; [p]). (A3) 

The third term is more complicated to derive since its expression is based on the exchange-correlation potential. 
Making use of the relation 



(c) 

the expression for h)J i(r becomes 



2.1 " \r-r>\ 2 J |r'-r"| 5p a (r) 



(A4) 



' dv"p ia {v") Sh ™ D ( r > f; . (A5) 



p(v) ^ " p(r) xc ^' — ^ ' p(r) J — ' <5p CT (r") " 
Regrouping all the terms, the expression for the exchange-correlation hole of the aNK functional can be written as 

NK — r „^ SD (r,r';[p£ f ]) 



^ c NK (r,r') = ^ c S °(r,r'; [p]) + _ £ (p 4(T (r)p(r') - p(r)p ta (r>) + p£ f (r) J dr"p la (r", 

+ / 15 (2/ rf -/ M )^(rK (r') + 4 SD (r I r'; [p - ^ CT ])(p(r) - p M (r)) - ^D( r , r ' ; [p])p( r ) + h^ D (r, r'; [p™ f ])p iCT (r)) . 

(A6) 

The result given in Eq. (1551) can be obtained from the above equation taking into account the xc-hole sum rule of the 
LSD functional, J dr'h^ B (r,r') — —1 (valid for any electron number). 
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